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ABSTRACT 

The branching ratios and CP asymmetries of B — > tttt decays measured in B factory experiments 
indicate a large ratio of color-suppressed (C) to color-allowed (T) tree contributions and a large 
relative phase between the tree and penguin amplitude. In this talk, I would like to report on a 
recent analysis to see whether the large C/T ratio can be explained within the QCD based model 
with i) a large contribution from the annihilation terms or with ii) large final-state-interaction 
phase between two different isospin amplitudes. We show that the current experimental data 
do not exclude either possibility but we may be able to distinguish these two effects in future 
measurements of direct CP asymmetry in B — > 7r°7r° decay. 



* Talk given at First Workshop On Theory, Phenomenology And Experiments In Heavy 
Flavor, Anacapri, Italy, May 29-31, 2006 



1 Introduction 



The large B — ► 7r°7r branching ratio measured at Babar and Belle indicates that the color- 
suppressed amplitude is large compared with predictions from the short-distance effective non- 
leptonic weak Hamiltonian. the ratio of the color-suppressed to color-allowed tree amplitude 
C/T is proportional to 02/01 with a\ = C\ + C%/N c and 02 = C2 + C\/N c being the effective 
Wilson coefficients for the charged and neutral current tree-level color-singlet operator 0\ and 
O2 respectively. At the renormalization scale [i ~ 1 , C\ ~ 1.10, C2 ~ —0.250 (at NLO) 
for non-leptonic B decay, 02 therefore becomes much smaller than a\. Thus in the absence 
of higher order QCD radiative corrections and non-factorisable contributions, B decays with 
neutral particles in the final state are suppressed compared to decays with charged particles 
in the final state and are called color-suppressed decays. This rule is known to be violated in 
D decays and also in B decays with charmed meson in the final state. Since we expect QCD 
Factorization (QCDF) 2 to work better for charmless two-body B decays with energetic light 
hadrons in the final state, the large C/T ratio extracted from the BaBar and Belle B — ► tttt 
measurements compared with theoretical predictions from factorization models is a surprise: 
The suppression of the B° —* tt + tt~ and the enhancement of the B° — > tt°tt° decays. Also 
the large measured CP asymmetry in B° — > tt + tt~ implies a large relative phase between the 
tree and penguin amplitude and suggests that strong final state interaction effect could play an 
important role in B — > tttt decays which are governed not only by weak interaction but also 
by strong interactions. Therefore an understanding of the strong interaction effects in these 
processes is crucial for extracting the weak phase and ultimately possible new physics which 
could make additional contributions to the one-loop penguin amplitude 3 . In this talk, I would 
like to report on a recent workjl] in which we investigate two possible mechanisms for the 
enhancement of C/T in B — * tttt decays: i) the higher order correction of QCDF and ii) the 
effect of final state interaction (FSI) and show that in QCDF annihilation terms could play 
an important role in the enhancement of C/T but alternatively, a large FSI rescattering phase 
could also enhance C/T in QCDF with moderate annihilation contributions. 



2 The decay amplitude from experiments 

With the recent measurements of the branching ratios for the B° — > 7r + 7r~, B° — ► 7r°7r° and 
B + — > 7r + 7r° as well as the mixing-induced S n + n - and direct C 7r + 7r - CP asymmetries in B° — > 
tt + tt~ decays, it is now pos sible to determine the ratio of penguin to tree amplitude and the 
relative phase between tree and penguin amplitude as done in many previous model-independent 
analysis [5]-[T8]. 

In the c-conventionjH], the B — > tttt amplitudes are 

A(B° -» vr+vr-) = Te i&T e^ + Pe iSp (1) 
V2k{B Q -> vrV) = Ce iSc e^ - Pe iSp (2) 
V2A(B + -> vr+vr ) = {Te i&T + Ce iSc )e i " / (3) 

where <5t,c,p are respectively the strong phase for the tree T, color-suppressed C and penguin 
P amplitude, 7 being the CP violating weak phase. To determine C/T and the relative phase 
between C and T, we have 5 measured observables |19| (all branching ratios quoted are CP- 
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averaged) 



c 


= -0.50 ±0.12 




= -0.37 ±0.10 


Br(vr + 7r") 


= (4.5 ±0.4) x 10 -6 


Br(7r°vr ) 


= (1.45 ± 0.29) x 10 


Br(vr + vr ) 


= (5.5 ±0.6) x 10" 6 



(4) 
(5) 
(6) 
(7) 
(8) 

B° mixing-induced and direct CP asymmetries 



where S„+ n - and C n + n ~ are respectively the B c 
in the expression for the time-dependent CP asymmetry of B 

A n+n -{t) = S n+7T - sm(AM B t) - C n+7T - cos(AM B t) 

In terms of T, C and P, we have: 



R S„ 



sin 2a + 2 sin(/3 — a) cos 5pt 



sin 2(3 



where 



R (' '- ! - =2 sin(a + (3) sin <5pt ( — 



/? - 1 - 2«>s(n + ?) ( ^) + 



The CP-averaged branching ratios are given by: 



2Br(7r°7T ) 1 



R+- 



C 



Br(vr+7r-) R 

2Br(7r + 7r°)r B o _ 1 
" Bi(tt+tt-)t b + ~ R 



—j ~ 2cos i s PT ~ &ct) cos 7 I — 



C 



I + 2cos<V / ( ^ + ^ X 



0) 

(10) 
(11) 
(12) 

(13) 
(14) 



To have an idea of the size of the relative phase 5pt, let us take as an approximation, a ~ 90° 
and neglect the (P/T) 2 term in Eg. 1)10(1 , we have 



tan 5pt ~ — CL 



(15) 



which gives 5pr — —36.5° more or less in agreement with the more precise value —41.3° for 
7 = 67° obtained in Table 1. This simple calculation shows that there is a large relative 
penguin-tree strong phase in B — » irir decays. Table 1 also shows the variation of P/T, R and 
5pt with 7 and CP asymmetries within the uncertainties in the measured CP asymmetries and 
the determination of 7. 

Figure ^ shows the allowed region for C /T and 5qt plotted with the central values of the 

CP asymmetries measurements. Note that the overlap of the solid (i?oo) an d the dashed (R^ ) 

bounds shift towards the larger C/T region as 7 becomes larger. R^ allows relatively small 

value of C/T while Rqo leads to a more strict constraint, C/T > 0.5. The overlap region is 
distributed in a large range of 5ct- Note also errors due to (S^^- , C 1T + 1T - ) are ± a few % in 
C/T and ±20° in 5ct- This analysis shows clearly that C/T is large compared to the leading 
order QCD corrections(C/T = 1/N C ). In the next section we present results of QCDF calculation 
for C/T. 
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(-0.62, -0.27) 


0.49 
-166° 
0.40 


0.29 
-147° 
0.69 


0.19 
-92.1° 
1.03 


0.32 
-43.4° 
1.35 


0.52 
-27.5° 
1.63 


0.72 
-20.5° 
1.82 


0.91 
-16.6° 
1.91 


(-0.50, -0.37) 


0.55 
-164° 
0.35 


0.38 
-149° 
0.62 


0.26 
-115° 
0.92 


0.29 
-66.6° 
1.21 


0.44 
-41.3° 
1.46 


0.62 
-29.8° 
1.63 


0.80 
-23.4° 
1.72 


(-0.38, -0.47) 


0.61 

-164° 
0.33 


0.45 
-150° 
0.58 


0.33 
-125° 
0.85 


0.31 
-86.9° 
1.12 


0.41 
-56.7° 
1.35 


0.56 
-40.4° 
1.51 


0.72 
-31.2° 
1.58 


(-0.38, -0.27) 


0.59 
-170° 
0.31 


0.40 
-162° 
0.55 


0.23 
-140° 
0.81 


0.17 
-79.1° 
1.07 


0.31 

-37.8° 
1.28 


0.49 
-24.1° 
1.44 


0.67 
-17.8° 
1.51 



Table 1: Determination of P/T (upper value), 8pt (middle value) and R (bottom value) using 
experimental results for 5' 7r + 7I — = (—0.50 ± 0.12) and C w + W - = (—0.37 ± 0.10) for given values 
of 7, 7 = (27° ~ 87°) . 



3 C/T in QCD Factorization 

Neglecting electroweak penguin contributions, T, C and P for B — > irir decays in QCD Factor- 
ization [2J HH! are given by: 

Te iST e^ oc X* u ( ai + h + a u A ) , (16) 
Ce iSc e^ oc A*(a 2 -6i-a£), (17) 
Pe iSp oc X*a c 4 (18) 

al = a\ + r x a p & + 26 4 + b 3 (19) 
r x = 2m^/(2m(,m 9 ) ~ 1.24 and m q = (m u -t-m^)/2 . 

a„ = v ub v: d ~ A\\p - i v ) (20) 
A c = V cb V c * d ~ -^IA 3 . (21) 

where ctj contain the factorisable term, the vertex correction, hard-scattering correction and 
penguin correction and 6j are annihilation contributions. Using the following parameters, 

H = 4.2 GeV, m q (2GeV) = 0.0037 GeV, A B = 0.35 GeV, 

|Au/A c | = \p + ir)\ =0.09, af = 0.1 (22) 

we find 

at = 1.02e io - 8 ° - Q.QU PH e^ H (23) 



3 




Figure 1: Allowed region for 5ct (x-axis) versus C/T (y-axis), from bounds for i?oo and R+-. 
Values for P/T and 5pt are from the central value of (£„.+„.- , C n + n -) = (—0.50,-0.37). The 
three solid lines represent Rqq = 0.64 — 0.14, 0.64, 0.64 + 0.14, the three dashed lines represent 

= 2.27 — 0.32, 2.27, 2.27 — 0.32. The overlap of solid and dashed bounds are the allowed 

region for C/T and 8ct- 



a 2 = 0.21e- i23 ° + 0mip H e i4>H (24) 
a\ + r x al = -0.097e i21 ° + O.OOlOpn e i0H (25) 
: -0.10e i7 ° + 0.0010p#e^ H (26) 



and 



bx = 0.027 + 0.063pV^ + 0.0085(p A e^ A ) 2 (27) 

6 3 = -0.0067 - 0.021/) A e^ A - 0.015(p A e^ A ) 2 (28) 

6 4 = -0.0019 - 0.0046p A e^ A - 0.0006( / 9 j4 e i<?iA ) 2 (29) 

We have made a complete analysis of C/T covering all the parameter space of p's and <^'s. In 
the limit of ph,a = . We find 

9± e i5cT Q = o.29e- i8 ' 5 ° (30) 
To 

Ph,a and 4>a,h represent the divergent end-point singularity terms for the hard scattering func- 
tion and the annihilation diagrams which introduce large theoretical uncertainties in the QCD 
calculation. Fig|5]shows 5qt versus C/T for — ir < 4>a,h < ^ with fixed pn and pa which would 
vary, say, in the ranges of < 1 ~ 2. We see that quite a large range of C/T and 8ct are 

allowed from QCD factorization, C/T up to 0.45 (0.55) for pa = 1(2). However the range of 
allowed values for C/T is still below the fitted values shown in Figure [I\£or 7 > 47° . We thus 
obtain a constraint 7 < 44° (52°) and 7 < 46° (56°), for pa = 1 and pa = 2 respectively. In 
|20| the problem of the small a 2 has already been recognized and a possible solution was pro- 
posed : by choosing the the largest value of the Gegenbauer moment of ir distribution amplitude 
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Figure 2: Scattered plot of the QCD factorization estimate for 5qt (x-axis) versus C/T (y-axis) 
including the end-point singularity effects. In the plot, we fix the p parameters as (phiPa) = 
(1, 1) (left) and = (1, 2) (middle) and vary the phases in the range of — tt < 4>a,h < tt (interval 
of 0.2 radian). The rest of the parameters are fixed (see text for details). The last figure (right) 
is obtained in the same manner with (j>h,Pa) = (1)1) but with different parameter set, the 
so-called scenario 2 of QCD factorization (see text for details). 




Figure 3: The left figure is the plot of Eq. (|34|). the result including the FSI phase, in the plane 
of S C t cB (x-axis) versus C cS /T cS (y-axis) by varying 5 2 o- versus C cS /T cS (y-axis) by varying 5 2 o. 
The number on the line indicates the value of 620 at each point. The bare value Cq/Tq = 0.29 
obtained from the default parameter sets of the QCD factorization with pn,A = is used. In 

the middle figure, we put together the left figure and the experimental bounds from and 

i?oo for the case of 7 = 57° (Fig.l upper-right). The right figure is obtained in the same way as 
the middle one but using Cq/Tq ~ 0.61, the result with the parameter set called scenario 2 in 
QCD factorization. 



(scenario S2). More recently, this model has been reanalyzed with NLO corrections to the hard 
spectator-scattering diagram and an enhancement by a factor of 2 for the hard-scattering 
correction is obtained: 

a 2 ~ 0.48e- il0 ° + 0.18p H e^ H (31) 

Co/Toe^o ~ 0.61e- i3 ° , p H>A = 0. (32) 

As shown in the scattered plot of Fig. 2 with p& = \ for scenario S2(right most), C/T ~ 1.1 can 
be achieved in this scenario if 5ct is very small. It appears that QCD factorisation can solve 
the large C/T puzzle. 



4 Does FSI Phase Make C/T large 



Recent studies show that FSI phase in B — > 7T7t can be large[22j and could effectively enhances 
C/T |23l I24j . since the -B — ► 7r°7r° decay can be induced by the charge exchange scattering 
process 7r + 7r~ — » 7r°7r° which can generate C from T. We examine here whether the FSI 
effect can enhance sufficiently the value of C/T of the QCD factorization without adjusting the 
incalculable parameters pn,A and <f>H,A- We now add the FSI phases to the QCDF amplitude 
and evaluate C/T and constrain the values of 820 = 82 — 8q and 5ct using the allowed region 
for C/T in Fig.l with 7 = 57°. In [23], a similar analysis is performed and a large 820 is found 
by a fit to the the central values of the experimental data. We have 

Tese^cff = [(2T - C )e iS ° + (T + C )e 4<52 ]/3 

Cefre^cff = [_(2To - C )e l5 ° + 2(T + C )e i52 ]/3 

P eS e iSp ee = P e i( - So+Sop \ (33) 

where Co, 7b, Po are C, T and P for = 

(£*L\ JSct s = (-2 + 2e^) + (l + 2e^)Co/r 

VTkJ C (2 + e «ao) + (-l + e i*«))C7o/ro 1 j 

(C/T) e ff indeed becomes larger as the FSI phase increases. The bare ratio Cq/Tq = 0.29 can be 
enhanced to (C cfr /T cff ) ~ 0.4 for 8 2 o ^ ±21°, where S C t cS — ±44°. 

For 7 = 57°, Co/To = 0.29, the allowed region from and Rqq overlap at £20 — —65°, 

at which (C/T) e g ~ 0.96 (middle figure). For 7 = 57°, Co/Tq = 0.61(scenario S2), the central 

value of ,Rqo) are reproduced by 820 — 40° where C e g/T e g ~ 0.8 and 8cT eS — 40°. To 

constrain 7, we need the FSI phase and input from QCDF. Then assuming ^20 30°, we obtain 
7 < 48° (55°) using the QCDF default values. With the parameters for scenario 2 shown above, 
we find that ^20 ^ 30° leads to 7 = (59±3)°(> 56°). Note that there may be a FSI contribution 
not only to the phase but to C/T itself [23], so the bounds obtained above may receive corrections 
from QCD or FSI effects. 



5 Conclusion 

Using the latest B — > irn data, we find that in B — > tttt decays C/T is large over a large range 
for Sct- C/T > 0.5 for 7 > 47°. We then show that QCDF with large pu.A leads to large values 
of C/T, especially when Sct is small. In a scenario of QCDF with ph,a = but with the FSI 
strong phase included, C/T and Sct are enhanced when the FSI phase £20 increases. A precise 
measurement of the B° — > 7r°7r° direct CP asymmetry Coo could distinguish these two models. 
In fact 

do _ C singer - 8 PT ) 1 ^ 
C^ T sin 8 pt Roo 

A small Sct(— 0) leads to a ratio of order unity with negative sign, Coo/C| ~ —(C/T)/R 00 , 

implying Coo = 0.57 for Sct = , close to the higher end of the current experimental value 
Coo = 0.281q39 while a large 8ct(~ =' =7r /2) would imply a strong dependence on Spt for 
the ratio C 00 /C+_ ~ ±(C/T)/(R 00 tanS PT ). 
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